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INTRODUCTION 

The assumption that the input source or calling popu- 
lation is a Poisson process is common to many operations 
research models such as waiting lines, inventory systems, 
and maintenance schedules. When this assximption cannot be 
supported, computational techniques assuming Poisson input 
are inapplicable and the scope of the operations, analyst ' s 
study of the system may be severely limited. 

Epstein^ gave a fairly comprehensive review of methods 
for testing the validity of the assumption that the under- 
lying distribution of life is exponential. Reinmuth^ sug- 
gested a "simple statistical test for detection of a Poisson 
input source" based on the fact that events generated by 
such a source are uniformly distributed over time, i.e., the 
probability that an event will occur within any given interval 
of time remains constant. This particular test was among 
those reviewed earlier by Epstein. Scheuer and Trueman^ dis- 
cussed several methods of testing including Reinmuth's test 


and also the modified Kolmogorov-Smirnov test described 
by Lilliefors^ for the exponential distribution with mean 
unknown . 

Any of these methods would suffice for many applica- 
tions provided sufficient data were available for testing 
purposes . Since valid data may be extremely limited for 
many practical situations^, a computer simulation study 
was carried out to compare the power of the Reinmuth or 
Z test with that of the Kolmogorov-Smirnov or KS test for 
the exponential distribution against a wide range of 
alternatives . 


SIMULATION PROCEDURE 

In order to carry out the study under conditions which 
would be applicable to a wide range of practical problems, 
two classes of alternative distributions were used. The 
first was generated by substituting uniformly distributed 
random nxjmbers for the cumulative distribution function of 
the two parameter Weibull and solving for the inverse. By 
varying the constant for the shape parameter, the nature of 
the resulting distribution of times between events could be 
varied widely for different cases. Thus, for a shape 
parameter of 3.3 the distribution is approximately normal. 


2 


I 



whereas for a shape parameter of 1.0, the distribution is 
exponential . 

A second class of alternatives was sought for which 
the nature of the distribution remained constant but the 
rate parameter changed gradually as a function of time. 

Such a class of distributions should be useful for repre- 
senting the daily or seasonal variations in the rates of 
customer arrivals at serving lines, etc. A cursory survey 
of the simulation literature was carried out but failed to 
suggest an algorithm for generation of data with a changing 
rate parameter. However, after some deliberation, it was 
noted that times between events could be generated for any 
distribution function, these values could be laid out in 
sequence on a time line starting at some appropriate point 
to give times of occurrence of successive events, these 
times could in turn be transformed by use of some appro- 
priate mathematical function, and finally a set of trans- 
formed times between events could be obtained by subtraction. 
For the present study, the exponential distribution was 
selected, the starting point was taken to be time equal to 
one, and the transformation was effected by simply raising 
the times of occurrence to some exponent in the range of 0.4 
to 2.0. Values for this exponent of less than unity gave 
decreasing times between successive events whereas values 
greater than unity gave increasing times. Times generated 
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in this manner were considered to represent a Poisson type 
process with a gradually changing rate parameter. 

A simulation case was considered to consist of all data 
generated for some specific value of either the shape param- 
eter of the Weibull distribution or the exponent used to 
generate a Poisson process with changing rate parameter. 

Each case consisted of 5000 replicate r\ins and provided data 
for a nxjmber of different sample sizes . To minimize the 
effects of sampling variations, the random number generator 
was set to the same initial value at the beginning of each 
simulation case and was not reset tmtil the next case. Dif- 
ferent initial values were used for the two classes of alterna 
tive distributions. 


ANALYTICAL PROCEDURES 

For the KS test, the test statistic given by Mann, Schafer, 
and Singpurwalla^ was calculated for each sample size of each 
simulation run and compared with the critical values given by 
Lilliefors . 

For the Z test, the test duration is normally some speci- 
fied period of time rather than some specified number of events . 
However, for this study, the duration of each portion of each 
simulation run was determined by the time at which the last 
event occurred plus a correction amounting to one half the 
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average time between events for that run. Z values were 
calculated in accordance with Reinmuth's procedure and 
compared with the critical values of the normal distribution. 


VALIDATION OF SIMULATION RESULTS 


Lilliefors gave simulation results for the power of the 
KS test when the alternative distribution was log normal. 

To obtain an overall validation of the simulation procedures 
used in this study, data for the log normal distribution was 
generated using the function suggested by Hahn and Shapiro . ^ 

The results for the KS test are in good agreement with those 
of Lilliefors as shown below. 

Sample Critical Probability of Rejecting Hypothesis* 

Size Level Lilliefors "Ihis Study 


10 

.01 

.023 

.037 

10 

.05 

.082 

.089 

20 

.01 

.046 

.046 

20 

.05 

.113 

.136 

50 

.01 

.085 

.096 

50 

.05 

.215 

.255 


*Both studies based on 1000 simulation runs. 

RESULTS 

Tables 1 and 2 give the results for the power of the KS 
test against alternative distributions consisting of the two 
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parameter Weibull with shape parameters ranging from 0.5 to 
4.0. Note that for a shape parameter of 1.0, the distribu- 
tion is actually exponential and the results are in excellent 
agreement with the expected values of 1 and 5 percent for 
sample sizes up to around 70. For larger sample sizes, the 
results exceed the expected values by amounts which appear 
to be significant. Although this could indicate a bias in 
the random niamber generator, it could also indicate that 
the Lilliefors^ approximation method for calculating critical 
values is slightly biased for the larger sample sizes. Since 
the difference between the observed and expected values was 
small, no attempt was made to determine the exact cause. 

Results for shape parameters other than 1.0 indicate 
that the KS test is quite powerful for this class of alterna- 
tive distributions with only seven instances for which the 
power of the test was less than 80 percent at the 5 percent 
level of significance. Note that for a shape factor of 3.3 
for which the Weibull distribution approximates a normal 
distribution, the power of the test was 94.3 percent for 
the 5 percent level of significance and a sample size of 
only 10. 

Since Reinmuth's Z test is dependent on a changing rate 
parameter rather than the nature of the distribution function, 
it would not be expected to reject the null hypothesis -for 
this class of alternative distributions. Values determined 



for the power of the Z test for a few simulation cases 
served to confirm this expectation and are not presented. 

Results for the KS test are given in Tables 3 and 4 
for the second class of alternative distributions. Note 
that when the exponent used for transformation is equal to 
one, no transformation takes place and the resulting distri- 
bution is exponential. Inspection indicates that the KS 
data for this case are, as expected, very similar to those 
given in Tables 1 and 2 for the Weibull distribution with 
shape parameter equal to one. 

Results for the Z test are given in Tables 5 and 6 for 
this same class of alternative distributions. Values for 
the power of the test when the exponent was one are somewhat 
higher than the expected values for the smaller sample sizes 
To determine if this anomaly was due to the approximation 
method of determining test duration used for this study, a 
number of additional cases were simulated using an exponent 
of one and varying the correction added to the time of occur 
rence of the last event. The results, not shown, were sensi 
tive to this variable and indicated that a correction of 
approximately 0.4 gave results in close agreement with the 
expected values for sample sizes of ten. 

Inspection of the results in Tables 3 through 6 for 
exponents other than one indicates that the Z test is con- 
sistently more powerful than the KS test against this class 



of alternatives with differences greater than 70 percentage 
points being determined for some cases . 


CONCLUSIONS 

For the class of alternative distributions represented 
by the two parameter Weibull, the KS test for the exponential 
distribution is quite powerful whereas the Z test is worth- 
less. On the other hand, for the class of alternatives 
with changing rate parameters, the Z test is consistently 
more powerful than the KS test with differences in excess of 
70 percentage points for some cases . It is therefore con- 
cluded that both tests should be used for applications for 
which detailed knowledge regarding the possible classes of 
alternative distributions is lacking. 
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TABLE 1. POWER OF THE KS TEST AT 5 PERCENT LEVEL OF 
SIGNIFICANCE WHEN THE ALTERNATIVE DISTRIBUTION WAS 
THE TWO PARAMETER WEIBULL 

Values in the body of the table represent probabilities (in percent) of 
rejecting the null hypothesis. 
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TABLE 2. POWER OF THE KS TEST AT 1 PERCENT LEVEL OF 
SIGNIFICANCE WHEN THE ALTERNATIVE DISTRIBUTION WAS 
THE TWO PARAMETER WEIBULL 

» 

Values in the body of the table represent probabilities (in percent) of 
rejecting the null hypothesis. 
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TAWLK 3. POWKR OF THE KS TEST AT 5 PERCENT LEVEL OF SIGNIFICANCE WHEN THE 
ALTERNATIVE DISTRIBUTION WAS A POISSON PROCESS WITH CHANGING RATE PARAMETER 

Values In the body of the table represent probabilities (in percent) of rejecting the null hypothesis 


0.4 I 0.5 


Exponent Used to Obtain a Changing Rate Process 


10 
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5.6 
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20 
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8.8 

6.8 

30 
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12.4 

B.O 

40 

31.9 

16.5 

m 

50 

41.3 

21.4 


60 

51.2 

27.0 

12.9 

70 


15.0 

60 

68.5 

37.1 

16.4 

90 

75.6 

I 

1 

100 
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5. 

0 

5.2 
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m 
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.5 
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.1 
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.±_ 

CM 

.7 
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00 
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TABLE 4. POWER OF THE KS TEST AT 1 PERCENT LEVEL OF SIGNIFICANCE WHEN THE 
ALTERNATIVE DISTRIBUTION WAS A POISSON PROCESS WITH CHANGING RATE PARAMETER 

Values in the body of the table represent probabilities (in percent) of rejecting the null hypothesis. 
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TABLE 6. POWER OF THE Z TEST AT 1 PERCENT LEVEL OF SIGNIFICANCE WHEN THE 
ALTERNATIVE DISTRIBUTION WAS A POISSON PROCESS WITH CHANGING RATE PARAMETER 

Values in the body of Che Cable represent probabillcles (In percent) of rejecting the null hypothesis. 
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